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Abstract 

We investigate the interaction of ultracold antihydrogen with a conducting surface. Our discus- 
sion focuses on the physical regime where the phenomenon of quantum reflection manifests. We 
calculate the reflection probability as function of incident atom energy. We find that, for ground 
state H atoms (with T < 10^^ K), the probability of refiection is ii ~ 1 — kb, where k is the momen- 
tum of the atom and b = 2174.0 a.u. is a constant determined solely by the van der Waals-Casimir 
tail of the atom-wall interaction. 

We show that quantum reflection, which suppresses the direct contact of ultra-cold atoms with 
the surface, allows for the possibility of conflnement and storage of cold antihydrogen atoms. We 
calculate the life-time of confinement as a function of antihydrogen energy. We develop a theory 
of H in a wave-guide and propose its application to fundamental measurements. In particular, 
for measurement of retardation corrections in the long-range component of the antiatom - wall 
potential. We demonstrate, for H falling in the gravitational field of Earth onto a conducting 
surface, the existence of quantized H states. We calculate that the lifetime of ultracold H in its 
lowest gravitational state and obtain r = {Mgb/2h)~^ ~ 0.1 s, where Mg is a gravitational force 
acting on the antiatom. We propose that measurement of this lifetime may provide a new test for 
the gravitational properties of antimatter. 

PACS numbers: 
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I. INTRODUCTION 



The phenomenon of quantum reflection in the ultra-cold atomic collisions has attracted 
lots of attention during the last decade. The propensity of an ultra cold atom, under the 
influence of an attractive atom-wall potential, to under^ total reflection has been predicted 
on quantum- mechanical grounds quite some time ago ill. However, the observation of this 

[in 

counter-intuitive effect has only recently been observed in the laboratory [2], p[. Because 
quantum reflection prevents atoms from reaching the wall, it suppresses inelastic atom- 
surface reactions and thus, it can be exploited as an important tool for manipulating cold 
atoms Another reason for interest in quantum reflection is that it occurs at large atom- 
wall distances where retardation effects are important. Retardation leads to the Casimir 

effect BBS 

, where the van der Waals atom-surface potential melds into a 1/2;^ power 
law potential at large distances. Quantum reflection is sensitive to the long-range component 



of the atom-surface potential 

QED flflQIn]. 



and may therefore allow new tests for the predictions of 



Recent success in the production of the cold antihydrogen (further referred to as H) 
Q, [3] has spurred renewed interest in this fleld. The ultimate goal of the antihydrogen 
project, initiated at CERN, is to enable accurate tests of CPT and gravitational properties 
of antimatter. Such a capability requires the availability of ultra-cold antiatoms. Quantum 
reflection may provide a new tool that hastens the realization of laboratory measurements. 
We will show that an interesting manifestation of quantum reflection, for antiatoms, is the 
existence of long-lived states of ultracold H in material cavities. We introduce a theory 
of wave-guides for anti-atoms and discuss their possible application in laboratory studies. 
Measurement of anti-atom life-times in a cavity may provide valuable information on the 
properties of the long range interaction between the antiatom and the surface. The mea- 
surement of the gravitational force acting on an antiatom poses as an intriguing possibility. 
An essential property of antiatom-surface interaction is annihilation of H on the surface. 
Anti-atoms that are not reflected are lost to annihilation and H-wal\ annihilation events can 
be detected to measure the quantum reflection probabihty. 

In Sec. II we discuss quantum reflection of H impinging on a perfectly conducting surface 
at normal incidence and present a calculation for the reflection probability. The latter is 
given as a function of incident energy, including the zero energy limit characterized by the 
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complex scattering length. In section III we study the probabihty for quantum reflection as 
function of the antiatom - wall distance. In Sec. IV we apply our theory and calculate the 
lifetimes of H enclosed in a material cavity. In Sec. V we discuss the passage of H through 
wave-guides. In Sec. VI we discuss the possibility for application of quantum reflection in 
measurements of the gravitational properties of an antiatom. Finally, in Sec. VII we study 
the loosely bound states of atoms conflned in the atom-wall potential. We show that such 
states strongly affect quantum reflection when absorption by the wall is weak. We discuss 
the physical conditions corresponding to the limiting cases of weak and strong absorption 
(annihilation) on the surface and show that, in the latter case, the (anti) atom- wall states 
are destroyed. 



II. QUANTUM REFLECTION AND ABSORPTION 

Quantum reflection from an attractive (with a sufficiently large derivative) potential is 
known to manifest 14| at low energies. In our discussion we are mainly interested in those 
cases where the probability of quantum reflection is close to 1. The qualitative parameter 
associated with the quantum motion of the particle with mass m and energy E in the 



potential V{z) is the local de Broglie wavelength Xsiz) = 2TTh/ -^ImiE — Vi^zy). One 
expects quantum behavior when: 

> 1. (1) 

For a homogeneous potential, —Cg/z^ with s > 2 in the limit of zero energies condition ([T]) 
is fulflUed for distances z such that 

Zc < z^Zf (2) 



= {2^2^Jsf'^'~^^ (3) 



Zf = l/V2mE. (4) 

Condition ([2]) deflnes the domain where the reflected wave is generated and is applicable 
if the energy of the incident atoms satisfles the condition 

E<^Ec = l/{2mzl). (5) 



A detailed discussion of the ultra-cold atom quantum reflection can be found in refs. [8|, Il5|, 



16|. The H - wall interaction is described by the potential, hereafter denoted V{z). It is 



not homogeneous and differs from that of the H - wall interaction. In particular it is purely 
attractive even at very short distances z < Zg — 1 a.u. and is strongly absorptive at the 
origin due to the likely-hood of annihilation. However at large H - wall separations z ^ Zs 
it is dominated by induced dipole-dipole terms in the atom-wall interaction, similar to the 



case of the hydrogen - wa 
wall is known for large z 



1 system. The potential for H impinging on a perfect conducting 
and will be hereafter denoted as Vcp{z); in our treatment we 
use Vcp{z) calculated in [l7|. At distances Zs <^ z <^ X^^ the potential has the van der Waals 
form 

Vcp{z) ~ -Cs/z^ (6) 

with C3 = ^ ad{iuj)duj = 0.25 a.u.^. Here ad{iuj) is the dynamic dipole polarizability of 
the (anti)hydrogen atom, expressed as a function of imaginary frequency iu, and A^^ is the 
effective wavelength that gives the main contribution to ad{ioj). For the distances 2; ^ 
retardation effects are important and the potential is given by 

Vcpiz) ^ -CJz^ (7) 

where C4 = J;:^^^^ = 73.62 a.u.^, ad(0) = 9/2 a.u.^ is the ground state static dipole 
polarizability of the antihydrogen, and a = 1/137.04 is the fine-structure constant. We note 
;hat the above expressions for C3 and C4 are valid in the limit of perfect conducting surface 
5], and can be specified for a realistic metal or dielectric surface at finite temperature |l8|. 
At short distances -2 ~ -2^ the interaction between H and the wall differs from that of H. In 
particular it includes an inelastic component, corresponding to the process of capture of p 
and e in the medium of the wall followed by subsequent annihilation. We take annihilation 
into account by imposing full absorption on the wall implemented in the manner described 
below. 

The Schrodinger equation that governs the H-wall scattering is 



V{z) - E 



^{z) = (8) 



2mdz'^ 

At distances Zs z the potential has form and the solution of Schrodinger 

equation is 

^z) ~ (^Hl'\p) + exp(2i5,)i7f V)) (9) 
p = 2^/2mC3/z (10) 



where H[^\p) and H[^\p) are the Hankel functions [19] of order 1, and Sg = 5i + i52 is a 
complex phase-shift produced by the short range part of the interaction. The imaginary 
part of 6s is due to annihilation of H on the surface. We used here the analytical form of the 
zero-energy solution in homogeneous potential —C^/z^ and neglected the energy of the 
incident H since it is small compared to the potential Vcp{z) at 2; ~ ^g. In our treatment 
we are interested in low energies where the above requirements are always satisfied. 
The condition of full absorption, which we apply, is 

S2 > 1. (11) 

This condition selects the "incoming wave" solution of the Schrodinger equation and sup- 
presses the "reflected wave" 

$(z) ~ v^ijj'^(p) , <^<A,. (12) 

One can check that in the limit of small z the above solution coincides with the semiclassical 
"incoming wave" 

^/^H[^\2^2MC3/z) ^ ^ exp(-z / p{x)dx) (13) 

^/p{z) J 

where p{z) = ^J^mC'^l z^ is the local classical momentum. 

The full absorption condition ^2 ^ 1 leads to insensitivity of solution (|T2l) in details of 
the short-range interaction since it is independent of the short-range phase shift bg- Thus, 
for strong annihilation on the surface the outgoing flux of li is determined solely by the 
asymptotic properties of the potential. 

We impose (1121) as a boundary condition in solution of the Schrodinger equation, which 
now does not depend on the short-range physics. The large z asymptotic form of such a 
solution is 

$(2: — s> 00) ~ exp(— ifc^;) — S'exp(ifc2;) (14) 



where k = y/2mE is the momentum of incident H atoms and S is the diagonal element of 
the 5— matrix that describes elastic scattering. S defines the flux reflected by the asymptotic 
potential and the reflection coefficient R = jS'p has been obtained by numerical integration 
of eq. ([8]) with eq. f|T2l) as a boundary condition. The resulting reflection coefficient is 
presented in Table [T] as a function of the incident energy. It is seen from the last column 
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that the reflection becomes considerable (> 50 %) at ^ 10~^^ a.u. (T ~ 2.7- 10~^ K ) and 
reaches 99 % at ^ 10~^^ a.u. (T ~ 2.7 ■ 10"^ K). The results underscore that, in spite of 
strong atom- wall attraction, slow (anti) atoms are reflected by and only a small fraction of 
the incoming flux reaches the wall. 

The reflection coefficient in the low energy limit is conveniently expressed in terms of the 
energy independent constant a = liuik^o^^, known as the scattering length. In the low 
energy limit the amplitude of the reflected wave 5* can be expanded in terms of k 

S{k) ^ 1 - 2ika. (15) 

Due to annihilation on the wall l^l < 1, the scattering length a acquires a negative imaginary 
part. The flux reflected in the region where the long range tail of the interaction potential 
dominates takes the form frequently used in the literature on quantum scattering 8| 

|5p ~ 1 -4A;| Ima| = 1 - A;6 (16) 

where we have introduced the constant b = 4|Ima|. The flux, absorbed (annihilated) on 
the wall is 

P{k) = 1 - \S\^ ~ kb. (17) 

The annihilation probability P becomes small when kb <^ 1, whereas the reflection proba- 
bility tends to 1. We have calculated the scattering length for the potential Vcp{z) with full 
absorption boundary condition, 

= -52.4-^543.5 a.u., (18) 
Reacp/lmacp = 0.09, (19) 
b = 2174.0 a.u.. (20) 

The annihilation probability P for potential Vc-p is expressed in terms of the scattering 
length and is tabulated in Table [T] along with the exact values. One notices from column 
4 that the scattering length approximation (eq. (ITTl) becomes valid for E < 10^^^ a.u. i.e. 
when kb <^ 1, but the validity of the non-perturbative approximation (given in column 3) 
extends to higher temperatures. 

It is interesting to compare this value for the scattering length with the corresponding 
value for the purely homogeneous potential —C^/z"^ with the full absorption at the origin. 
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log(-E/a.'u.) 


P 


1 — exp(— fc6) 


kb 


R 


-9 


0.95 


0.99 


4.16 


0.05 


-10 


0.69 


0.74 


1.32 


0.31 


-11 


0.33 


0.34 


0.42 


0.67 


-12 


0.12 


0.13 


0.13 


0.88 


-13 


0.04 


0.04 


0.04 


0.96 


-14 


0.013 


0.013 


0.013 


0.987 


-15 


0.0042 


0.0042 


0.0042 


0.9958 


-16 


0.0013 


0.0013 


0.0013 


0.9987 


-17 


0.00042 


0.00042 


0.00042 


0.99958 


-18 


0.00013 


0.00013 


0.00013 


0.99987 



TABLE I: The annihilation (P) and reflection (R) probabilities for the ultra-cold antihydrogen 
impinging on the material wall. Column 3: non-perturbative approximation to P. Column 4: 
scattering length approximation to P. 



In the latter case 21 1 



af' = exp(-Z7r/(s - 2)) 



s-2 



2/(s-2) 



r((.-3)/(.-2))/r((. -!)/(.- 2)) (21) 



and for s = 4 (the Casimir correction) the scattering length is purely imaginary and given 
by 

04 = -i\/2mCi (22) 

which for H results in 

a^{H) = -i519.9 a.u.. (23) 

We note that the imaginary part of the scattering length acp is rather close to the value 
obtained for a purely homogeneous —6*4/^^ potential and suggest that the main contribution 
to the penetration probability arises from the asymptotic Casimir tail ([7j). The nonzero real 
part of the scattering length is the contribution from distances where the potential changes 
from the van der Waals ([6]) to the Casimir limit (JTj). In fact, one can see from ( |2T]) that 
the real part of the scattering length for the purely homogeneous potential with s = 4 (with 
absorptive core) is exactly zero, while for s < 4 it is negative. 
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We investigate the influence of tlie inner van der Waals part of the potential in the zero 
energy limit. We model the exact potential Vcp{z) by the analytically solvable potential [8| 



with / = C4/C3. Such a potential has correct asymptotic behavior in the limit of big 
[z ^ /) and small {z <^ /) distances. The zero energy wave-function for such a potential is 



$™(^) = v/i(7+7) (iff)(2v/2mC3(l/z + C3/C4)) - e^'''H^'\2^2mC3{l/z + C^/C,))^ 

(25) 

To obtain the scattering length we examine this expression in the limit z ^ I to first 
order in terms proportional to l/z 

<!>m{z) ~ z{l + l/{2z)) (Hi'\ai + i-)) - e'^'^H^iiil + ^))) (26) 

where ^ = 2C3A/2m/C4. Employing a Taylor expansion and collecting terms proportional 
to z, the wave function becomes 

Now we use the fact that the asymptotic expression of the wave-function can be written 
in terms of the scattering length as 

^„,{z > /) ~ 1 - z/ara (28) 
Comparing the last two equations we extract the scattering length from eq. 

The full absorption boundary condition (^2 ^ 1, e^^^" — > 0) leads to, in analogy with the 
results of f[T^ . cancellation of all terms that include the solution H^\^) 



The ratio Ream/Imam is given by 



Reajlma^ = ^ + ^^(0 + ^AiOMO + YliOYiiO)) (31) 



where Ji{z) and 1^1(2;) are Bessel functions 19| of order 1. For H we obtain 



= -69.8 - i505.6 a.u., (32) 
Re Om/ Imam = 0.14, (33) 
e = 3.53 (34) 

The ratio Re a^/ Imam contains important information on the scale parameter ^ that, in 
turn, determines the transition between the van der Waals and Casimir limits. In the limit of 
large ^ the model potential Vm becomes —C4^/z'^ whereas the ratio becomes Ream/ Ima^ — 
1/^ as it should for the purely homogeneous potential V = —C^/z^. According to this 
treatment, we expect that for the exact potential Vcp, the smaller the ratio Re a^/ Imam, 
the more important the contribution from the Casimir tail to the scattering amplitude 

III. PROBABILITY OF QUANTUM REFLECTION AS A FUNCTION OF DIS- 
TANCE 

We now study the relative importance of different regions in which the van der Waals- 
Casimir potential dominates and contributes to the value of the reflected wave. We represent 



the wave-function in the following form 



22| 



^{^) = -I /I A ( B{z) exp(i [ p(x)dx) — exp(— z f p{x)dx)\ (35) 

where A{z) and B{z) are arbitrary functions that are obtained when we substitute $(2) into 
the Schrodinger equation. Here p{z) = ^j2m{E — Vc'p(2;)) is the local classical momentum, 
and Zq is an arbitrary distance. Identification of a solution of the Schrodinger equation with 
the incoming or outgoing wave is unambiguous only in case when it can be represented in the 
semiclassical form exp(±i J^^p{x)dx). If the the WKB approximation is valid in the entire 
range of antiatom-wall distances, the solution everywhere has the form of an incoming wave 
and no reflection occurs. However, for certain regions ([2]) of antiatom-wall separation, the 
WKB approximation fails (the so-called badlands). The solution of the Schrodinger equation 
differs from the semiclassical one and leads to the appearance of a reflected wave. Hence 



the function B{z) in (135|) can be interpreted as a function that "converts" the semi-classical 
solution into the exact quantal solution and thus contains information on the reflected wave 
amplitude generated at each separation distanee z. 



We need to put one more condition on A 



'z) and B{z) in order to define them uniquely. 



Following the phase function method [22, l23| we require that 

^'{z) = iA{z)p^^'^{z) (B{z)exp{i f p{x)dx) + exp{-i f p{x)dx)\ . (36) 

V J ZO J ZO J 

As we have already mentioned, the function Bi^z^ is an amplitude for the reflected wave. The 
relationship between the asymptotic value B{oo) and the S-matrix can be easily established 

pz — *oo 

S = B{oo) exp(2i( / p{x)dx — kz)). (37) 

J ZO 

Substituting eq. (l35ll and eq. ( l36il into the Schrodinger equation we obtain the first order 
nonlinear differential equation for B{x): 

B'{z) = [B^{z)exp{2i [ p{x)dx) - exp{-2i [ p{x)dx)] . (38) 

\ J ZO J zn J 



As long as in the limit z ^ Zs the wave-function should have the form given by (1121) and 
( fT3l) one should choose the initial condition as B{zs) = 0, and in equation ( l38l) one should 

put Zo = Zs- 

The function can be interpreted as that "portion" of the reflected wave that is 

generated in the domain between at zq and z. (Such an interpretation is unambiguous 
when z belongs to the range where the WKB approximation is valid, and the functions 
exp(±2 J^^p{x)dx) can be identified with the reflected or transmitted waves). Obviously 
|i?(oo)p gives the "full" reflection probability for the given energy. In Fig.l we plot the 
function |i?(z)p for two energies E = 10~^^ a.u. and E = 10~^° a.u. 

As z —>■ the WKB approximation becomes valid. Using it (and thus reflection from these 
parts of potential exponentially decreases) we found that the contribution to the amplitude 
of the reflected wave from the distances z < 100 a.u. is indeed small. 

In addition, at ultra- low energies (practically for E <^ 10~^° a.u.) the reflected wave is 
generated within a wide range of distances. In that case we cannot define an unambiguous 
"reflection distance". The plot shows that for the energy E = 10~^^ a.u. the reflection 
probability is 0.88 and this value is in harmony with the calculated value presented in Table 
m We note that 75% of the reflected wave amplitude is generated in the domain from 500 a.u. 

10 




FIG. 1: Contribution of different antiatom-wall distances to the reflection probability. |S(2:)|^ 
expresses the reflection probability accumulated in the interval between zq and z. 

to 5000 a.u., the region where the Casimir tail of the potential dominates. At intermediate 
distances, from 100 a.u. up to 500 a.u., the contribution is 22%. However, with increasing 
energy the reflecting domain is more localized and is shifted to shorter distances, whereas 
the reflection probability is diminished. According to Fig.l, for E = 10^^^ a.u. the reflection 
probability is 0.31 and 92% of the reflected wave is generated in the domain from 100 a.u. 
up to 1000 a.u. In this case the reflection probability reaches its full value within a well 
defined interval Az around Zr, and allows us to define Zr, the "reflection distance". At higher 
energies, reflection occurs at the reflection distance ( Az <^ Zr ) and the reflection probability 
becomes exponentially small. The reflection coefficient as well as reflection distance z^ can 



be estimated using "complex turning point method" |8|, l2J, |25 |. 

As discussed above, the real part of the scattering length is more sensitive to the details 
of the potential at intermediate distances. More generally, one can expect that the real part 
of the phase shift S{E) (defined through the equality S = exp{2iS{E))) provides additional 
information on the antiatom-surface interaction. In Fig. [2] we plot Ke6{E) for the exact 
interaction Vcp{z), as well as for two model interactions Vi = —Ci/z'^ and V2 = —C^/lz'^ + 
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FIG. 2: HeS/k as a function of energy for exact interaction Vcp (solid line) and for model inter- 
actions Vi (dashed line) and V2 (dotted line). 

We found that the exact phase shift is negative for E > 10~^ a.u. and positive for 
E < 10^^ a.u. The phase shift produced by the homogeneous potential Vi = —C^/z^ is 
always negative and tends to in the limit E 0. Note that in the limit of zero energy 
Re6{E ^0)/k = -Reacp. 

In conclusion, we emphasize that annihilation of H on the surface allows a description 
that depends only on the nature of the van- der Waals-Casimir potential. In particular 
all inelastic processes that occur at characteristic distances of few a.u. (including thermal 



heating of H atoms by phonon exchange 



26 



271]) do not affect the reflected wave. 



In the case of strong absorption the phase and amplitude of the reflected wave are in- 
dependent inn details of the short-range interaction. This allows one to glean important 
information on parameters that characterize the van der Waals-Casimir atom-wall potential. 
In the limit of a perfectly conducting surface, they are solely determined by the dynamical 
dipole polarizability of H. Our treatment can be generalized to the case of finite conductiv- 
ity or non-zero temperatures by introducing corrections to the Casimir-Polder interaction 
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constants In particular, the simple estimation of H quantum reflection constant b in 
case of dielectric surface can be obtained using expression fl22|) with the modifled value of 
C4 coefficient: 

6^4W2mC4^^r(e) 
V e + 1 

where e is permittivity of the dielectric material. The correction function F (e) calculated in 
[g] is approximately equal to 0.77 in a wide range of e. The estimated in such a way constant 
b for polyethylene (e ~ 2.3) is 6 ~ 1138.9 a.u., which is approximately twice smaller than 
in case of perfect conducting surface. As a consequence the characteristic energy, for which 
quantum reflection from polyethylene surface starts to be signiflcant [kb 1, E = 1/ (2m6^)) 
is four times greater, than in case of perfect conducting surface. 

In the following section we consider the conflnement of antihydrogen between material 
walls. In particular we show how the properties of quantum reflection are related to the 
shifts and widths of the energy levels of H atoms enclosed in cavities with conducting walls. 



IV. ULTRACOLD H BETWEEN TWO CONDUCTING WALLS 

The partial reflection of ultra-cold H from the material wall enables the existence of 
quasi- stationary (decaying) states of H in-between two walls. In so far we will be interested 
in the one-dimensional wave-function of such states. We expect that the distance L between 
the conducting walls is large compared to \acp\ so that the atom interacts with each of the 
walls independently. The boundary condition at the left wall {z = 0) is given by f[T^ . To 
obtain the boundary condition at the right wall one should replace in equation ([? ]) z by 
L — z. Far from the walls (|acp| <^ z ^ L — \acp\) the wave function is 

$(2) ~ sm{kz + 5cp) (39) 

where 6cp is the phase-shift produced by interaction with the wall. Using the symmetry of 
the problem one obtains the quantization condition 

kL + 2Scp = vrn. (40) 

We consider the low energy domain where that the scattering length approximation is 
valid, so that 6cp = —kacp- Thus we get for k 
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k = nn/{L-2acp), (41) 

and so the energies of the box-hke states are quantized. The eigenvalues are modified by 
quantum refiection according to 



" - = 2miL - 2acp)^ " + " ^' '"^"^^ ^^'^ 

The widths of the states are given by 

^" = ''2;;^ = -IT 

where ei"^ = 7r^n^/(2mL^) are the energy levels unperturbed by the long range H - wall 
interaction. 

The level shifts and widths are determined by the value of complex scattering length 
[b = 4\Ima\, c.f. eq. [TSll that characterizes quantum refiection. The lifetimes increase with 
the wall separation as and decrease with the excitation quantum number as n^. These 
expressions are obtained using the assumption that the scattering length approximation f[T5|) 
is valid, and restricts its validity for e„ < 10^^^ a.u. 

The solutions ^n{z), corresponding to the complex energy levels Sn, are the eigen- 
functions of a non-self-adjoint Hamiltonian. In fact, they are decaying quasi-bound states 
with energy — ir„/2. This implies that such states obey the bi-orthogonality condition 

/■oo 

/ <^niz)^kiz)dz = Snk (44) 

Jo 

and differs from the "standard" expression by the absence of complex conjugation of the 
(bra| function. 

The expression for the width (1431) coincides with the simple formula for the annihila- 
tion rate of the particles moving freely with the velocity v = ^J^e'^^jM between the walls 
separated by the distance L 

r = P{v)uj, (45) 

where P[v) is the wall penetration probability and uo = v / L is the frequency of "hits" on 
the wall. Substituting the probability of annihilation on the wall P{v) from (|T71) we again 
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obtain the expression given in fl43l) . i.e. T = 2e^^^j^. At high energies {e > 10~^ a.u.) the 
annihilation probabihty P{v) ^ 1, and we get T = j^. The hfetime of such a "fast" particle 
is just its time of flight between the walls. 

As an example, we take L = 10 fim which is much greater than the length scale b = 0.115 
lim, that characterize the Casimir force. At this wall separation distance, the ground state 
energy is Eq = 7.5 10"^"^ a.u. the width is Fq = 1.7 10~^^ a.u. and corresponds to the lifetime 
0.014 s. We compare this value to the time of flight between the walls t = L/v = 0.0005 
s. Thus H in this state bounces about 30 times before it annihilates. Since the life time 
is proportional to L^, and the passage time to L^, the number of the bounces between the 
walls grows linearly with L. Let us mention here that the concept of quasi-stationary state 
itself is meaningful only in case its life-time is much greater than the corresponding time of 
flight between the walls. 



V. ANTIHYDROGEN IN A WAVE-GUIDE 



A laboratory demonstration of quantum reflection states could realized by introducing 
H atoms through a slit between parallel conducting walls separated by distance L. (The 
similar principle was used to observe quantum motion of neutrons in the gravitational field 



of Earth 



28|, 



291 ]. Ignoring, for the moment, the influence of gravitation on H, we align the 



conducting walls parallel to the gravitational field. The Schrodinger equation, that governs 
the motion of antiatoms inside the wave-guide with transverse dimension z and horizontal 
dimension x, is 



+ V{z) + V{L-z)-E 



^{z,x) 



0. 



(46) 



2mdx'^ 2mdz^ 

We express the two-dimensional wave-function "^{z, x) as a series of products of normal- 
ized transverse wave- functions ^n{z) given in eq. (139!) and longitudinal plane- wave functions 
exp{ipnx): 



^(2;, x) = Y^ Cnexp{ipnx)<l>n{z) (47) 

n 

where C„ are the amplitudes of corresponding states, dependent on the properties of 
the flux entering the wave-guide, and p„ is the horizontal momentum of the state with the 
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transverse energy Sn — iTnfi, so that 

pl/2m + - iTj2 = E. (48) 

Let us mention that the eigen-functions $„(2;) of not-self-adjoint Hamiltonian do not 
form a basis in the Hilbert space [sQ]. So far the expansion f HTI) is an approximation. We 
will show however that under certain conditions of the wave-guide experiment few quasi- 
stationary states from (1471) would be give exhaustive contribution to the measured flux, 
which verifies our approach. 

If the characteristic spatial dimension of the incoming flux Hq is much larger than the 
separation L and the distribution of the transverse velocity in the incoming flux exceeds 
the value 1/L, the first few excited states are uniformly distributed. Because the transverse 
energy — ir„/2 is complex, the horizontal momentum Pn, for a given state, is also complex 

(c.f. HHD. 

The horizontal momentum at small transversal energy (£„ <^ E) is given by 



„ ^ - V^^isn + tTj2) =p- '^"^^^ (49) 



where p = y2mE and v = p/m. We can thus write for the wave- function 'if{z, x) 

'^{z,x) = exp{ipx) Cn exp(— r„a:/(2f )) exp —ienx/v^n{z). (50) 

n 

Using the definition t = x/v we relate the decay time of the transverse states 
$„(z) exp(— r„t/2) to the disappearance rate (along the x-axis) of the flux density inside 
the wave-guide. 

The integrated flux density F = Fq J |^(z, d)\^dz (where Fq is a normalization constant) 
gives the number of counts of H in the detector at the exit {x = d) of the wave-guide 



F = F^Y. |Cn|'($n|$n) exp(-r„r^'«^'') (51) 

n 

(V -\-V \'T-P"'^^ 

+ FoJ2 C:C,{^n\^k) exp(^(£fc - e„)r^'«-) exp(- ^ " ^ ) 

where we have introduced the passage time r^"** = d/v. 
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We denote the second term in this expression as an "interference" term since its presence 
is due to non-orthogonahty of the decaying states, (^^l^^fc) 7^ 5nk- The appearance of such 
a term is not surprising. Indeed, decaying states are not stationary and they don't have well 
defined energies and so they can be thought of as a superposition of truly stationary states. 
Therefore transitions between states with frequencies ujnk = — occur. The transitions 
manifest in the appearance of the oscillating exponents exp(ic<j„fcrP"'''^) in the expression 
for the measured flux flSTl) . However, observation of such interference requires rather good 
resolution in the horizontal velocities of the initial flux. The uncertainty in passage time 
^rj-pass _ ^pass^^i^ -g order of and it follows that the limit for horizontal velocity 
resolution is 

- < T^'^^'Unk. (52) 
V 

For a broader distribution in the horizontal velocity, the interference terms cancel due to 
fast oscillating exponents and the expression for the flux takes the form 

F = Fo5^exp(-r„r^'"-). (53) 

n 

The time of flight r^*^*** may be chosen in such a way that only a few transverse states 
with the smallest r„ = n^Ti contribute to the flux of H in the wave-guide. This time of 
flight selects the states with small transversal energy and thus our expression for the 
horizontal momentum (H9l) is justified. Let us also note that such a choice of time of flight 
ensures that contribution of quasi-stationary states to the flux at the exit of the wave-guide 
is much greater than the contribution of all the "non- resonant" terms, which were neglected 
in the approximative expansion ( 1471) of the total wave-function [31]. The typical length d 
for a flux of H with the horizontal velocity f = 10 m/s should be about 10 cm to ensure 
that only the ground transverse state passes through the slit L = 10 /im. Changing the 
slit separation leads to an enhanced contribution from the excited transversal states, due to 
the dependence (H3|) of the decaying rates r„ on L. The measured function F{L) provides 
information on the parameter h. In particular, if only the ground state passes through the 
guide, one gets the following (Ci and Fq independent) expression for h 

d\nF{L) \ mL^ 

dL ) 37r2rP"'^^' ^ ' 
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We stress that the quantity b above is very sensitive to the Casimir force since quantum 
reflection in the guide selects only those antiatoms that have interacted with the surface 
through the van der Waals - Casimir potential (all others are annihilated). 



VI. GRAVITATIONAL EFFECTS 

We have shown that quantum reflection takes place for very slow anti-atoms. This fact 
could be exploited to study the gravitational properties of anti-atoms. In this section we 
discuss the possibility for measuring the gravitational force acting on antihydrogen. We dis- 
tinguish between the gravitational mass, which we refer to as M and inertial mass, hereafter 
denoted by m. We consider an H atom bouncing above a conducting surface in the gravita- 
tional field of Earth. Confinement is achieved by quantum scattering by the van der Waals- 
Casimir potential from below, and by the gravitational field from above. The characteristic 
length and energy scales are 



The characteristic scale of the gravitational energy oi H is e = 2.211 10~^^^ a.u., whereas 
the corresponding length scale is Iq = 5.871 fim and is much greater than the length scale, 
(6 = 0.115 fim ), of the Casimir force. Thus, the gravitational field does not significantly 
affect the atom-wall interaction. The energy levels and the corresponding wave-functions 
are 



En = e\n, (57) 
^n{z) = CAi(^//o-A„) (58) 

where Ai(x) is the Airy function |32|!] and the eigenvalues \n are found by enforcing 
the condition Ai(— A„) = 0. The modification in the values for the eigenvalues A„ due to 
quantum reflection are obtained by matching the wave-function of H reflected from the wall 
(which at the distance \acp\ ^ z <^ Iq has the asymptotic form $(2;) ~ 1 — z/acp) to 
the gravitational wave-function Ai(2;//o — A), where A is a modified eigenvalue. Taking into 
account that in the matching region 2;//o ^ 1, we obtain the following equation for A 

18 



= -acp/lo. (59) 

Ai'(-A) 

Because acp/h <^ 1, A can be expressed as a perturbation series in terms of the gravita- 
tional eigenvalues A^ 

An = A„ + 5\n (60) 
with 5\n ^ 1- This leads an expression for the modified eigenvalues A„ 

An = A„ + acp//o- (61) 
From it, we obtain modified energy levels and the corresponding widths 



i?„ = e(A„ + ^), (62) 

Tn = e J- (63) 

The widths of the gravitational states fl63|) are proportional to the combination s/Iq. 
Using (1551) and (156!) we find that this ratio is equal to the gravitational force s/Iq = Mg so 
that 

r„ = ^-Mg. (64) 

The value of the ground state energy of H in the gravitational field of Earth turns out 
to be El = 5.17 10~^^ a.u. The corresponding life time is 

r = — - ~ 0.1s. (65) 

Mgb ^ ^ 

We note factorization of gravitational effects (appearing in the above formula via fac- 
tor Mg) with the effect of quantum reflection, appearing through the constant b. Such a 
factorization is a consequence of the small ratio in the characteristic scale b/lo ~ 0.02. 

We find that the widths of the gravitational states ( l63l) are independent of the energy (for 
the energies En < 10^^^ a.u.) to the first order in the small ratio 6/A„. This is understood 
using the following simple argument. The frequency of the atom bouncing above the surface 
in the gravitational field is a; ~ 1/ \/E, whereas the probability (|T71) of annihilation on the 
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wall is P ~ v-E'- Combining the two, we obtain an energy independent expression for the 
width, r = uoP. That is, all H atoms with < 10^^^ a.u. bouncing on the surface survive 
for 0.1 s. 

The measurement of this lifetime can be achieved by monitoring the spatial decay of the 
flux of antiatoms with the velocity v parallel to the surface. Such a flux, as function of the 
distance d and integrated over vertical dimension z is 

F{d) = Fo^M-\^fn = ^M'^r^n- (66) 

If one shifts the position d of the detector and measures the spatial decay F{d), one can 
infer the gravitational force Mg acting on the antiatom since 

d\nF(d)hv . 



VII. SURFACE STATES FOR ATOMS AND ANTIATOMS 



Our treatment of quantum reflection of antiatoms makes use of the fact that the antiatom- 
wall interaction is strongly absorptive at short distance, and is mathematically expressed by 
absorption condition (fTTIl . In this section we study the case where absorption on the surface 
is not necessarily strong. Such a situation may correspond to the partial loss of ordinary 
atoms that undergo inelastic interaction at the surface, or to reflection of antiatoms from 
partially transparent evanescent-wave atomic mirror. 

Below we outline a resonance-type dependence of the atom loss probability on the proper- 
ties of the short range interaction in case of the weak surface absorption. We show that such 
resonance behavior stems from the presence of narrow near-threshold atom-surface states 
33|. 

For the purpose of illustration, and in order to allow an analytic description, we replace 
the exact van der Waals-Casimir potential by the model potential (|2^ . We do not impose 
the full absorption boundary condition, instead the complex phase-shift 6s, produced by the 
short-range interaction is used. 

For the model potential Vm, the analytic expression for the scattering length (129|) simplifies 
substantially since ^ = IC'^^lmj C/^ ^ 1. and using the asymptotic representation of the 
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Hankel functions h[^'^\^) ~ y^2/{7r^) exp{±i{^ — Sn/A) we obtain 

a^ = -^ecot(e-5i-z52-37r/4). (68) 

We explore the dependence of on the real part of the short-range phase-shift 6i. We 
find that am is an oscillating function of 6i. The smaller the imaginary part 62, the more 
prominent are the oscillations. The maxima occur when 

Sl^' = ^ - 37r/4 - vrn, n = 1, 2, ... (69) 

This expression is simply the condition for the appearance of a new state in the potential 
well of an attractive van der Waals- Casimir potential and a short-range core. 

The existence of the narrow near-threshold state leads to a probability loss P = Ak\lmam\ 
that varies rapidly with variations in the phase Si. At its maximum, the probability loss 
(for the above mentioned case when ^ ^ 1 and 62 <^ 1) is 

Pma. = 2klj- (70) 
O2 

and the width of such a maximum is equal to 62- 

This example illustrates the influence of near-threshold singularities on the S-matrix and 
scattering observables. To find the position of the S-matrix poles in the complex /c-plane, 
we express the S-matrix in the vicinity of its pole 

S=i±|. (71) 

z — k 

Taking into account that for small k <^ \z\ we have S = 1 — 2ika, we get for the near- 
threshold pole position 

z = i/a. (72) 

Using a Taylor expansion of the scattering length a^ fl68l) in the vicinity of the resonance 
Sl^"^, the position of the near-threshold pole can be written as 

(73) 

where A61 = Si — Sl'^^ is the deviation of the short-range phase-shift from the resonance 
value. 
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The maximum in the loss probabihty occurs when the pole z of the S-matrix crosses the 
real axis of the complex /c-plane. In the absence of absorption, the pole lies on the imaginary 
axis of the complex fc-plane and when it crosses the real axis the scattering length approaches 
positive infinity. In the presence of absorption the pole has a real negative component. The 
smaller the shift from the imaginary axis, the more pronounced the effect. 

As an illustration, we consider the values 5i and 62 obtained within the model potential 

^ri \ } t^exp(-«(^)exp((2;o-2;)/p) \i z>zq 

Vs{z) = { (74) 
V exp if z < zq 

Here the argument of the complex exponent relates the real and imaginary part of the 
complex potential and p is the diffuseness radius. This model allows analytic solutions and 
facilitates our understanding of inelastic processes nearby and inside the wall. 

For a strong complex exponential potential Vs where | Im p^2M'U\ ^ 1, the scattering 



length is 



34| 



;2o + P (27 + 2 ln(pV2m|f/|) - (75) 



where 7 ^ 0.577 is the Euler constant. 

We first consider the case where the diffuse radius p of the exponential potential is much 
smaller than the local de Broglie wavelength X^iz^ = 1/ \/2mC-il z\ at the distance zq- Here 
the short range potential Vg varies rapidly when compared to variations of the de Broglie 
wavelength. At distances p ^{z — z^ ^ Ab(2;o) the wave-function is 

^, = 1- (76) 

Os — Zq 

Matching it with wave-function ( 125|) and taking into account the asymptotic properties 
of the Hankel function with a large argument 

eo = 2 v/2mC3(l/zo + C3/C4)) » 1 (77) 

we obtain 



h^i^- 37r/4 (78) 
52 = <^p/Ab(^o) = VP\l'^MC^/zl (79) 
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where we have kept only the leading terms in p/Xb{zq). 

Thus, the limit p/Xb{zq) implies vanishing absorption (^2 — > 0) and corresponds to 
the case where the reflecting wall is placed at position zq. 

The weak absorption limit ^2 ^ 1 can be realized by an optical evanescent wave that 
generates an effective repulsive potential. The detailed study of this case is beyond the 
scope of our paper; for the purpose of our qualitative analysis it is only important that the 
imaginary part of the phase-shift 62 turns to be small when p/\b{zq) ^1. At the same 
time the phase-shift 5i is well approximated by expression (!78|) . 

In the case p/Xb{zq) > 1, the short range potential varies slightly with respect to the 
scale of change of and one can use the semiclassical expression for the phase-shift 



Thus a shallow but extended imaginary potential can produce a large imaginary phase- 
shift 62 ^ 1- A a necessary condition for the validity of the strong absorption boundary 
condition is that the short-range interaction varies slowly in comparison with the variation 
in the de Broglie wavelength Xb{zs) at Zg, the distance where the short-range interaction is 
important. For H the characteristic de Broglie wavelength at Zg — 1 a.u. is Xb{zs) ~ 0.05 
a.u. and the scale p of the short-range inelastic interactions is about 1 a.u. Well within the 
range where the strong absorption boundary condition is valid. 



where n = {^q — ^)/7i is the number of atom-surface bound states. 

The scattering length exhibits resonance behavior when the position of the "reflecting 
wall" is chosen in such a way that a new state appears at threshold, i.e. C,q{zo) — ^ — > irnQ. In 
the vicinity of zq it is useful to represent the scattering length ( ISTI) in terms of a Breit-Wigner 
form 




(80) 



We substitute (JZHD into (EHD to get 




(81) 



Im a, 
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(82) 



2 (e - ^0 - v™o)2 + 



Rea- 



2 (e-eo-vrno)2 + <5i' 



(83) 
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FIG. 3: Imaginary part of the scattering length as a function of the position of the "reflecting 
wall". 
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FIG. 4: Real part of the scattering length as a function of the position of the "reflecting wall". 
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In Fig. [3] we plot the imaginary part of function of zq for two values of the 

imaginary phase-shift 62 = 0.1 and 62 = 1.0. We vary zq from 10 a.u. to 15 a.u. Fig. H] 
illustrates the behavior of the real part of the scattering length for the two values of 62 given 
above. We find a resonance-like increase in atom loss when the short range potential is 
centered around Zq = 12.1 a.u. This resonance is very prominent for weak absorption when 
62 = 0.1, and it is practically washed out in case of stronger absorption 62 = 1.0. 

In the hypothetical case where the short range part of the (anti)atom-wall interaction 
can be varied in a controlled way (e.g. by reflection from the optical evanescent wave) one 
could accurately measure the properties of the (anti) atom- wall potential by careful tuning 
of resonances and fitting the experimental data to eq. fl82p . In this way, one could measure 
the values of the C3 and C4 coefficients. 

In conclusion, we find that the main difference in quantum reflection in cases of strong 
and weak absorption by the material wall is that in the latter case, narrow atom-surface 
states can be formed. The position of these states depends on the short range interaction. 
The presence of a near-threshold state strongly enhances the probability of atom loss, and 
might be used for resonant spectroscopy (e.g. with the use of optical evanescent wave). In 
case of strong absorption, the near-threshold states are completely destroyed by annihilation. 

VIII. CONCLUSIONS 

We presented a theoretical account of cold antiatoms interacting with a material wall. 
We focused on the antihydrogen - wall system where the latter is represented by an ideally 
conducting surface. We have shown that, as antiatoms approach the surface in an attractive 
potential they are not necessarily annihilated. If sufficiently cold, they may be totally 
reflected. Reflection occurs in the region where the long range tail of the van der Waals - 
Casimir interaction dominates. We calculated the reflection probability as function of the 
incident energy of antiatoms. We analyzed the behavior of the reflected wave as a function 
of distance from the surface and introduced the notion of a "reflection distance". We have 
also shown that quantum reflection allows the existence of long-lived metastable states of 
H in the conducting cavity. 

We pointed out that scattering of ultra-cold H in an attractive van der Waals - Casimir 
potential is sensitive to retardation in the antiatom-wall interaction. We explored the possi- 



25 



bility of performing experimental measurements with H by introducing it into a wave-guide. 
The rate of the H annihilation on the wall of the guide could provide important informa- 
tion on retardation effects in matter- antimatter systems. In the strong absorption limit, 
such measurements are not "contaminated" by effects produced on the wall and which are 
characterized by poorly known short-range parameters. However, they could be employed 
to study the asymptotic properties of H- wall interaction. 

We also proposed measurement of effects due to gravity acting on the antihydrogen. In 
particular we found that the life-time of H, in the combined gravitational and the mirror 
image potential of a conducting plate, to be on the order of 0.1 s. The established dependence 
of that lifetime, r = 2h/ (Mgb), on the gravitational force allows the possibility for measuring 
the gravitational properties of antimatter. 

Finally, we outlined the role of the near-threshold atom-surface states the where absorp- 
tion on the surface is small. For antiatoms, weak absorption could arise with the introduction 
of an evanescent wave that prevents antiatoms from direct contact with the wall. We showed 
that a resonant increase in the absorption rate can be observed if the atom- wall interaction 
is "tuned" to form an atom- wall state at the threshold. We note the possibility of exploiting 
resonant absorption for studies of atom-surface interactions. 
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